Hamiltonian structure of the Sabra shell model of turbulence: 
exact calculation of an anomalous scaling exponent 
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We show that the Sabra shell model of turbulence, which was introduced recently, displays a 
Hamiltonian structure for given values of the parameters. As a consequence we compute exactly a 
one-parameter family of anomalous scaling exponents associated with 4th order correlation functions. 
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The field of turbulence and turbulent statistics cannot 
pride itself on a large number of exact results. One of 
the best known exact result is the "fourth-fifth law" of 
Kolmogorov, which pertains to the third order moment 
of the longitudinal velocity differences, fixing its scaling 
exponent £3 to unity Q . This follows from the conserva- 
tion of energy, a quadratic invariant, by the Euler part 
of the Navier-Stokes equations. This famous result reap- 
peared also in simplified toy models of turbulence, like 
shell models, since the conservation of a quadratic in- 
variant was built into their definition [Q. Nevertheless, 
the high degree of simplification involved in the shells 
model did not lead so far to additional exact results. 

In this Letter we report a discovery of a Hamiltonian 
structure of the Sabra shell model of turbulence that was 
introduced recently The Hamiltonian structure ex- 
ists for given values of the parameters of the model, and 
it does not coincide with the usual quadratic invariants; 
in fact it is cubic in the velocities. There is a Hamil- 
tonian density TL n which is local in shell-space, and its 
existence implies a flux of a local conserved density. The 
consequence of the constancy of this flux in shell-space 
fixes the value of a fourth order scaling exponent which 
is "anomalous" in the sense that it must involve the exis- 
tence of a renormalization scale. This result is both exact 
and nontrivial, constituting a kind of "boundary condi- 
tion" on theories of anomalous scaling in such models. 

The Sabra model of turbulence, like all shell models of 
turbulence, describes truncated fluid mechanics in wave- 
number space in which we keep N + 1 "shells" . Denoting 
the nth wave-number component of the velocity field by 
v n , < n < N, the model may be written as 

v n = -i(ak n+1 v n+ iv n+2 - bk n v n ^iv* n+l (1) 
-cfc n _iu*_ 1 u n _ 2 ) - vk 2 n v n + /„ , 

where the coefficients a, 6, and c are real, v is the "vis- 
cosity" and f n is a forcing term, usually limited to acting 
on the first shells only. Equation ([!]) is equivalent to the 
Sabra model suggested in || after the change of variables: 
v n = u n for n even and v n = — tt* for n odd. Like in all 
shell models one builds in the conservation of "energy" 



that a + b + c = 0. The actual values of the wave- numbers 
k n are determined by fco and the "level spacing" param- 
eter A, k n = ko\ n . By rescaling we can always choose 
a = 1, leaving us with only two free parameters in this 
model, A and e = —b. The inviscid and unforced part of 
our model then reads 



v n = -ikn+i \y n +\v n +2 + eA~ 
+(1 - e)A~ 2 v„_ 2 <_i] 

Like the popular GOY model || 
conserves the so-called "helicity" 



N 



Vn-lV n+ i 



(2) 



the Sabra model also 



(3) 



For the special choice of e equal to the golden mean 
e 9 = (V5-l)/2, Eq.(D reads 



n+l 



A 



Vn+lVn+2 + y w n-l«„+l 



(4) 



V n -2V 



n-1 



By a direct calculation one can verify that in this case 
the model exhibits a third integral of motion which we 
denote as W: 



W = 



N-l 

E 

rn— 



] e g k Q [v n -xV n Vn+i + c.c] 



(5) 



The main statement of this Letter is that the Sabra 
model (||) with e chosen at the golden mean realizes a 
Hamiltonian structure in the sense that the equations of 
motion can be written in a canonical form 



dU{a v a*} 
i — 

da* 



(G) 



N 

n=0 



in the inviscid limit v — > by requiring 



with H{a,j, a*} being the Hamiltonian function of the set 
of pairs of canonical variables {oj , a*} representing all 
the degrees of freedom. The Hamiltonian H. is precisely 
the cubic invariant W, expressed in terms of canonical 
variables a n , a* , which are related to the velocity as fol- 
lows: 



f 



a n = v n /t n g ' 2 . (7) 

According to (||) and (0) the Hamiltonian has the form: 
ft = E!=i ftm where 

ftm = eg^o(A^") m [a m+ ia m a*„_ 1 + c.c] . (8) 

The statement is proven by verification: compute dji ac- 
cording to the canonical equations (0) with the Hamilto- 
nian TL: 

m 



a n = -i-Q-f = -*fcoe(A^")"[A % /i^a„ +2 a n+ i (9) 



+a 



,+i<2n-i 



:a n _ia ri _ 2 J 



Using the relationship ^ this equation coincides with 
Eq. (@), Q.E.D. 

One should note that in the majority of physical sys- 
tems in the conservative limit the Hamiltonian, if it ex- 
ists, coincides with the energy. Examples abound: elec- 
trodynamics, hydrodynamics, acoustics, waves on the 
surface of fluids, spin waves etc, see for example (J]). 
There are no regular procedures to find the canonical 
variables in terms of the "natural" variables of a par- 
ticular physical system. This relationship may be very 
nonobvious, like the Clebsch representation of the veloc- 
ity field in hydrodynamics. There are also examples of 
physical problems with a Hamiltonian structure in which 
the Hamiltonian does not coincide with the energy, like 
induced scattering of electrons on ions in strongly non- 
isothermal plasmas [^). The present example belongs to 
this class. 

The existence of a Hamiltonian has significant conse- 
quences for the solutions of nonlinear problems. In this 
Letter we demonstrate one such consequence, which is 
the ability to determine exactly an anomalous exponent. 
The locality of Tl n in the space of shell numbers means 
that there exist a current Ji{k n ) over shells of the mean 
value of the "density" (H n ). Here the symbol (...) stands 
for an appropriate ensemble average. In order to de- 
rive an expression for J 4 {k n ) compute the time deriva- 
tive of Tin (||) with the help of the canonical equations 
(||). The right-hand-side of this equation is a combina- 
tion of 4th-order correlation functions. The conservation 
of the Hamiltonian means that X^ft™ = 0- This is pos- 
sible only if (ft„) can be written as a difference of cur- 
rents depending on two neighboring shells (substituting 
the divergence in the continuous limit): 



(H n ) = J 4 (fcn-i) - J^K) 



We find 



J 4 (fcn) 



3 fc 2 (— F 4 (k n ) 
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(10) 



(ii) 



where Fi{k n ) is the following combination of 4th order 
velocity correlation functions : 



Fi(k n ) = Im[A 2 (w*_ 1 u„w n+2 w„ +3 ) (12) 

+£g{v* n -2 v n~lV n +lV n+ 2) 

In the steady state (ft n ) = 0, implying that J 4 {k n ) must 
be independent of k n : Ji{k n ) — Ji(k n -i). According to 
( 1 lh we must demand 



F 4 (k n ) ex ) 



(13) 



or, looking for a scaling solution F 4 (k n ) oc fc n ^ 4 , we com- 
pute 



a = i ogx ©= 2 + 1 o gA (^±I) 



(14) 



This is an exact solution for a one-parameter family of 
scaling exponents as a function of the level spacing A. 

We note that the anomalous exponent C, 4 > 2. The 
Holder inequalities imply that the scaling exponent £4 
characterizing any of the correlation functions constitut- 
ing F 4 {k n ) is smaller or equal to 4/3. Denote the scale at 
which the Hamiltonian invariant is pumped into the sys- 
tem by k n . If this scale is chosen to coincide with the first 
shells of the model, the anomalous exponent found here 
will always be subleading. This means that the leading 
scaling contributions exactly cancel in the combination of 
correlation functions making up F 4 (k n ) due to the con- 
servation of the Hamiltonian. On the other hand, if k n is 
chosen in the bulk of the inertial interval, and if this has 
an effect on the region of smaller k < k n , then in that 
region the contribution of £4 will be leading. This and 
other consequences of the existence of the Hamiltonian 
structure will be investigated further in the near future. 
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